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t^; Abstract 

' We review ongoing research related to the asymptotic dynamics of isotropic universes in theories 



with higher derivatives, especially near the initial singularity. We treat two major cases, that 



00 ■ is universes in vacuum, and also those filled with radiation using the method of asymptotic 

O ! 

, splittings of vector fields. Our solutions have the form of formal Fuchsian-type series for the 

CSJ , basic unknowns and are valid near collapse or big rip regimes asymptotically. 

o 

m 

Isotropic universes in vacuum or those filled with radiation have always been primary examples 

in efforts to model the effects of higher derivative terms near extreme conditions and cosmological 

^ ■ singularities, cf. Refs. [1], [2], [3]. With the advent of newer and more powerful mathematical 

5-j ■ 

methods to study the asymptotic structure of dynamical systems of Refs. [U [5], there is now the 
prospect of understanding the possible effects near spacetime singularities in a deeper and more 
complete and transparent way. 

Radiation-filled universes in gravity theories with higher derivatives have been analyzed in this 
way in Refs. [GjjT], in an effort to understand their dynamics near the initial cosmological singularity. 
The general higher order action is 



S = 1; [ Aotald/ig, (1) 
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where /^totai is the lagrangian density of the general quadratic gravity theory given in the form 

Aotal = + -^matter, with 

C{R) = R + ^R^ + 7Ric2 + (JRiem^, (2) 
where /3,7, (5 are constants. The field equations derived are as follows: 

(3) 
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2RR^''' - ^R^g^"" - 2{gf'Pg'"' - gf'''gP'')S/pS/aR 



where we have set = 2(3/3 + ^ + 6). This naturally splits into 00- and n-components (i = 1, 2, 3), 
but only the 00-component of ([3]) will be used below. Using the standard FRW metric, the field 
equation ^ leads to our basic cosmological equation in the form 



k + d 



a2 



'd d dd^ d'^ /c^"' 



(4) 



where is a constant defined by the constraint 



^ = ^, (from V^T^o = 0). (5) 

and k is the (constant) curvature normalized to take the three values 0, +1 or —1 for the complete, 
simply connected, flat, closed or open space sections. 

Setting X = a, y = d and z = d, Eq. dU can be written as an autonomous dynamical system of 
the form 

x = ffc,RAD(x), x = (x,t/,z), (6) 

where the curvature-radiation vector field ffc,RAD : M'^ — > : {x, y, z) fk,KAD{x, y, z) that defines 
the system is given by the form 

. , / - k"^^ 3y^ yz y k ky\ 

f.,RAD(x,y,z) = [x,y, + ^ + ^- -- ^- ^ + ^J- 

The curvature-radiation field ffc,RAD combines the effects of curvature and radiation and describes 
completely the dynamical evolution of any radiation-filled FRW universe in higher order gravity. 
Following the asymptotic methods of Ref. [1], we conclude that the only acceptable asymptotic 
splitting of the vector field f ^"^I^aq is of the form 

f _ f (0) If (sub) /ON 

t fc.RAD - t R AD + t fe,RAD ' i^i 



with with dominant part 

2x'^ 2y X 



faAD(x) = (y,^,;^ + ^-^), (9) 
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and subdominant part 

(sub) / C"^ ~ k'^S, y k ky 

ffc,RAD(x) = 1^0, 0, 2^^2y - ^ - ^ + ^ 

We then look for the possible asymptotic solutions, asymptotic forms of integral curves of the 
curvature-radiation field ffc,RAD; that is we search for the dominant balances determined by the 
dominant part f^^j^^D gi'^sn by Eq. Q. For this purpose, we substitute in the dominant system 
(x, y, z)(t) = f^'^I^AD forms x(t) = atP = {OfP, -qf^, pf) and solve the resulting nonlinear algebraic 
system to determine the dominant balance (a, p) as an exact, scale invariant solution. This leads 
to the unique curvature-radiation balance B^^rad € x Q'^, with 



^fc,RAD = (a,p) = 




where is a real, arbitrary constant. In particular, this means that the vector field f ^ had ^ 
scale invariant system. 

After calculating the spectrum of the Kovalevskaya matrix which is of the form 

spec(/Cfe,RAD) = {-1,0,3/2}, (12) 

in order to obtain the the /C-exponents of the dominant system, we proceed with the construction 
of the final solution by substituting the forms 

oo CXD oo 

x{t) = Y,ciiti^K y{t) = ^C2iti--2, z{t) = Y,C3iti-K (13) 

1=0 i=0 i=0 

where cio = 0,C2o = ^/2,C3o = —0/4:, and we are led to various recursion relations that determine 
the unknowns cij,C2i,C3j term by term. Further algebraic manipulations lead to the final series 
representation of the solution in the form: 

Our series ()14p has three arbitrary constants, 9, C13 and a third one corresponding to the arbitrary 
position of the singularity (taken here to be zero without loss of generality) , and is therefore a local 
expansion of the general solution around the initial singularity. Since the leading order coefficients 
are real, we conclude that there is an open set of initial conditions for which the general solution 
blows up at the finite time (initial) singularity at t = 0. This proves the stability of our solutions 
in the neighborhood of the singularity. 

The case of vacuum universes with higher derivatives is currently under investigation, cf. Ref. 
[8], both for the fiat and the curved subcases. The dynamics of the various families of models 
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depends on the asymptotic modes of the curvature-vacuum vector field 

ffc,vAc(2;, y, z) = (y^^- + - ~ lii ~ i^' ~^^^) ' ^^^^ 

and this leads to very interesting, possible dynamical regimes asymptotically. In the flat-vacuum 
case, we can show that there is a general solution of the field equations with the scale factor 
having the leading order behaviour t^/-^. In the curved cases, the situation becomes more involved 
and there can be various general and particular solution of great interest. However, it appears 
that the 'radiation-vacuum' solution mentioned above is a global attractor of the possible curved 
asymptotics. These results will be reported in detail elsewhere. 
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